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improved perturbative expansion in terms of unstable excitations. Our effec- 
tive propagators are derived from spectral functions with a constant width. 
These spectral functions describe screening and damping of gluons (photons) 
as well as "Brownian" motion of quarks (electrons). BRST-invariance allows 
us to reduce the number of independent width parameters to three. These are 
determined in a self-consistent way from the one-loop self energy and polar- 
ization tensor in the infrared limit thus rendering this limit finite. All spectral 
width parameters are found to be proportional to gT. We reproduce the well 
known expression for the electric "Debye" -screening mass. The transverse 
(magnetic) gluons (photons) are found to interact only at nonzero momen- 
tum or energy, at least to leading order. As a consequence their spectral 
function acquires a width only away from the infrared limit. Finally, plasmon 
modes are determined and found to be strongly damped. 



Work supported by GSI. GSI-Prcprint 95-4 (1995), subm. to Nucl.Phys. B 



1 Introduction 



It is well known that naive perturbation theory for gauge theories at finite tem- 
perature does not work. This is due to severe infrared divergences related to the 
approximation of massless gauge bosons having an infinite correlation length. The 
breakdown of perturbation theory can be understood qualitatively if one recalls that 
the average distance between "particles" in a heat bath is of order 1/T, the inverse 
temperature. On the other hand, due to the uncertainty relation, each "particle" 
has a quantum mechanical correlation length 1/M, where M is its mass. This im- 
plies that the individual wave functions substantially overlap as soon as T > M. For 
massless fields this condition is met already at infinitesimally small temperature. 

As an example providing a more concrete condition consider the gluon self energy 
at cij = 0,p ^ 0. The Z-loop contribution can be estimated as p. 145] (Lorentz 
indices are suppressed) 



where / > 2, (7 is the coupling constant and M is a fictitious mass of the gluons 
serving as an infrared cutoff. Here, one can see the infrared divergences to become 
more and more severe the higher the loop order is. 

There exists at least a partial solution to the infrared problem. This is the 
method of "hard thermal loops" . The essence of this method is the introduction 
of thermal screening masses. For electric gluons one finds a mass of the order of 
gT. This suggests that the loop terms in ([|) are suppressed by factors of g^^'^ as 
compared to the leading terms. Unfortunately, one finds no mass of this order for the 
magnetic or transverse gluons. This seems to imply that QCD at high temperatures 
finally is highly nonperturbative 0, |^ . 

With the present work we are presenting a slightly modified approach which may 
overcome these problems and possibly will allow for an expansion in the coupling 
constant. The main difference to conventional approaches is that we do not rely on 
quasi-particles, i.e., particles that have a definite mass-shell only slightly perturbed 
by the interactions. There are several arguments against the quasi-particle concept 
at finite temperature: 

1. A system at finite temperature does not exhibit the Poincare symmetry of 
the vacuum state. What remains instead is a semi-direct product of the four- 
dimensional translation group T(4) and the three-dimensional rotation group 




(1) 



spectrum 0, and only these representations may be used in a perturbative 
expansion 0. 



2. A self-consistent calculation of the spectral width of massive fermions at 
nonzero temperature T exhibits a nonanalytical behavior in T as T ^ 
1^. This is consistent with the picture of Poincare symmetry restoration at 
T = as a "thermal phase transition". It also tells us, that the basis for a 
perturbative expansion at T = is different from the basis at T > where we 
have to use a continuous mass spectrum [|], |]. 



3. Consequences of the KMS-condition 10| on quasi-particles have been inves- 



tigated by Narnhofer, Requardt and Thirring [|TT|. Their theorem essentially 
states that quasi-particles can be used as an approximation only if one consid- 
ers length and time scales larger that 1/T, i.e., the regime where the presence 
of thermal excitations is still negligible. 

4. Probably the most severe objection against the use of quasi-particles arises 
in the long-time limit, i.e., for times larger than f» 1/T. For such long times 
the assumption of a quasi-particle with an infinite life-time no longer holds 
because the presence of thermal excitations necessarily leads to scattering or 
absorption. Attempts to calculate the damping rates in a perturbative manner 
from the imaginary part of the self energy or polarization on-shell must fail 
if the mass of the quasi-particle is smaller than T. Indeed, calculations of 
damping rates indicate that one has to determine such rates in a self-consistent 
manner, although one frequently makes the assumption that the damping rate 



is much smaller than the mass of the quasi-particle [12, 13 



It is this breakdown of the quasi-particle picture for temperatures higher than the 
relevant masses on which we base our procedure. We are taking this feature into 
account establishing simple ansatze for spectral functions with a finite width. This 
will be the content of section 2. In section 3 we calculate the self energies and 
polarizations to one-loop order and determine from them the width parameters in 
a self-consistent manner. We show these results as well as those for plasmon modes 
and their damping rates for QCD, and in section 4 also for QED. Section 5 is devoted 
to a summary and conclusions. 



2 Spectral functions and imaginary time Green's 
functions 



The safest method for the determination of Green's functions uses dispersion rela- 
tions. They enforce the proper analyticity in the complex energy plane by removing 
unphysical poles from the physical Riemann sheet [0 . The central quantity of this 
procedure is a spectral function A{E,p) which must be known for real E. For 
arbitrary complex E, the retarded and advanced Green's functions read 

/OO 1 
JE'A(e;p)^-^^ . (2) 

This automatically implies that the spectral function is proportional to the imagi- 
nary part of the retarded respectively the advanced Green's function: 

A{E,p) = Tilm(^^'^(E,p)) = {S\E,p) - S'^iE.p)) . (3) 
71 zvri 

For perturbative calculations of the partition function one usually needs causal 
Green's functions in euclidean time respectively their Fourier transform S{pQ,p) 
with the Matsubara frequencies fl^ 



, yj" = 2ra ttT (bosons) . . 

Po={ : T f ■ \ ^ri = 0,±l,... . (4) 

y„ = [2n + Ij TTl (termionsj 

These functions can be obtained from the retarded or advanced Green's functions 
by analytic continuation in a unique manner [IIH IIH], using the relations 



S''{E,p) = ±S{po^e~iE,p), e-^0+ 

S^{E,p) = ±S{po^-e-iE,p). (5) 

With these relations we can now define S{po,p) via dispersion integral: 

r°° 1 

S{Po,p) = T dE A{E, p) —— . (6) 

J-oo E - ipo 

The upper sign in (j^) and (^) is valid for fermions, the lower one for bosons; they ac- 
count for the signs when switching to euclidean space-time and euclidean 7-matrices. 

We are interested in spectral functions with a finite constant width. For bosons, 
we abbreviate the spectral width as 7^ for the moment, and use the spectral function 



1 



AniUp \E — ijjp — i'-fB E — Up + i'jB 
^ 1 1 \ 

E + ujp + i-fB E + ujp - i-fB J 

where Up denotes the energy-momentum dispersion relation. We will neglect masses 

and simply choose uj"^ = p"^. The spectral function is normalized according to 

the commutation relations for bosonic fields: 

/oo 
dE' E'A{E',p) = 1 . (8) 
■oo 

The corresponding euclidean Green's function can now be obtained using eqns. (||) 
and ®: 

(bo I +7b) +P 

In the limit of vanishing width, 7^ — > 0, one recovers the usual perturbative Matsub- 
ara Green's function 1/(pq + p^). A remarkable feature is that the Green's function 
(H) is real, even with 7^ > 0. 

Similarly, we proceed for fermions, where we abbreviate the spectral width as 7^? 
to distinguish it from the Dirac 7-matrices (70,7)- As spectral function we take 

7^ 70(^2 +^2 + 72 ) + i2E7p 



Af{E,p) 



'joujp + iyp 'joUp + i-yp 



Aniujp \E — Up — ijF E — Up + i'-yp 

-foUJp - hp ^ 7o^p-i7P j j^^Q^ 



E + Up + ijF E + uop- hp 
with euclidean 7-matrices, see appendix A. This spectral function is normalized such 
that 

/oo 
dE'Tr[7o^(E',p)]=4 . (11) 
-00 

The corresponding Matsubara Green's function with euclidean 7-matrices and cj^ = 
p^ simply reads 

Q( ^ . 7osign(po) (bol +7f) +7P 

Sf{Po,P) = -i ^ — 2 ■ ^2 

(bol + iFr + p^ 

Obviously, the only modification as compared to the free case is a shift of the absolute 
value of the Matsubara frequencies by the spectral width parameter jp. 

For illustration let us consider the spatial and temporal behavior of a retarded 
real-time propagator as defined in (^). With a boson spectral function according to 
(0) we get 

/oo 



oc e(t)exp(-7Bt)^^ , (13) 

whereas the 70 component of the fermion propagator with spectral function ([TD|) 
becomes 

/oo 
dE exp(-iEt)Tr[7o^^(^,p)] 
-00 

oc 6(i(:) exp(— 7j7t) cos(c<jpt) . (14) 

This shows that 7 describes the damping of the corresponding wave-function. In 
order to investigate the spatial behavior we consider the time average of (pT^): 

/ rft5|(t,p) = -i5B(0,p) = ^— ^. (15) 

Obviously, for these zero frequency bosonic modes one cannot distinguish between 
a mass and a spectral width. This implies that a spectral width not only describes 
damping but also screening. 

Another interpretation of these spectral functions, which is easier to accept in the 
fermionic case, is the one of Brownian motion [p!7| , eq. (6.33)]. In this case, (p!3| ) resp. 
(p!4D represent the decaying probability amplitude to find a particle with momentum 
p in the same state after time t, when it is subject to thermal fluctuations of the 
surrounding medium. 



3 The partition function of hot QCD 

We now consider the partition function of QCD represented as a functional integral, 

Z = J VA^V^pVijVcVc exp{-SQCD) • (16) 



For convenience we consider a family of 0(3)-covariant gauges [|T8|. As we are 



interested in high temperatures we neglect current quark masses. The action reads 



^QCD 



+ i(A^A;(x))2 + K^e{x) {d,c\x) + gr'^A^ixYix))] (17) 

Here, /,m are flavor indices, are color indices and the field strength tensor is 
defined as 

F;, = d,Al-d,Al + gr'^AlAl . (18) 
The gauge fixing involves a "vector" A^ defined as 

A — /\fl J—fi\ /-\Q\ 



with gauge parameters ^ and A. 
are defined as follows: 



The usual perturbative 2-point Green's functions 



{i'i,a{x)ipm,i3{y))o = ^irr^S^p / jT^Z^^FQip) exp(ip(a; - y)) 



(2vr) 

Sfo{p) = {hoPo + hp)~' (20) 



for the quarks, 



mx)At{y))o = S'^' exp{ip{x - y)) 



(27r)4 



Sto{p) 
Ssoip) 
Siioip) 



Po + P 



pI + 



for the gluons and 

(c"(i)c'(y))„ = i"'|^S,„(p)exp(ip(x-»)) 

for the ghost fields. The symbol (. . .)o denotes functional integration with measure 
exp(— (So), taking into account only the bilinear part of the action. The integration 
over d^p^ is defined in appendix A. Recall, that for the corresponding objects in 
operator formalism time ordering is implied. 

In ( pTD we have introduced four symmetric, 0(3)-covariant tensors, 

pLl ^ . ,PiPi. 
' "m* p2 "l'' 



Among these, , and V^^ are projectors with 



V + V + P"^' = I4 . (24) 

The square of the traceless tensor V^^ is the projector onto the "longitudinal" 
subspace: 

= P^ + P^i . (25) 

Our aim now is to provide propagators having finite spectral widths. These spectral 
widths are assumed to be constant and will be determined self-consistently from 
the lowest order self energies and polarizations at zero momentum and energy. This 
means that we separate the unperturbed part of the action from the interaction part 
in the following way: 

Sqcd = So + Smt 

= 5o + Sy + Sint ~ S^ 
= So + Sint 

So = So + S^. (26) 

Perturbation theory is now done with respect to Sint, and So is the action of gen- 
eralized free fields, which already contain the full two-point correlations present in 
the theory. In particular, these generalized free fields may have a continuous mass 
spectrum, and therefore avoid the pitfalls of the quasi-particle picture we outlined 
in the introduction. This method has been introduced by Licht [0, its application 
to systems at finite temperature is discussed in refs. P, ||]. 

As this procedure constitutes a certain approximation scheme we have to worry 
about gauge dependence. The appropriate tool in this context are BRST-Ward 
identities. After gauge fixing the partition function of QCD has a global symmetry 
with Grassmann character, the so-called BRST-symmetry pO| : 

SipiAA = ic"(a;)T^,/3V^i,/3(a;) 
^d'ix) = -^r'^c'ixyix) 



This symmetry allows to derive the following Ward identity: 



{6[gA;{x)-c\y)]) = 

= {A';{x)AX{y) + d,c\x)c\y) + g Al{x)c\x)<^ {y)) 

(28) 

At leading order in an expansion in powers of the coupling constant this identity 
implies the following relations among the perturbative propagators, which also hold 
for the effective propagators of the generalized free fields: 

>^PoSe{,p) + ^\p\Sl2{.P) = PoSghip) 

\PoSl2{p) + ^\p\SLiip) = \p\Sgh{p) (29) 

These identities are fulfilled by the expressions in (^O]) - (^2]) as one can see by 
insertion. 

In order to simplify calculations we choose A = y/^ in the following. In use of the 
simple expressions we derived in section 2 for propagators with a constant spectral 
width, we now provide the following improved propagators for QCD: 

Sf{p) = (^lo sign(po)(bo| + If) + hP^ = 7m'^f,/.(p) 
1 

Srip) 



Se{,p) 
Sli{p) 



i\Po\ + Irf +P'^ 
1 1 

~a\Po\+iEY + \p^ 
1 



{\p^\ + ieY + \p^ 

Sl2{p) = 



Sghip) 



v^(bo|+7£;)' + |p^ 



(30) 



As these propagators describe generalized free fields [rU|, they fulfill the tree-level 
BRST identities (E^) by construction. We consider this a crucial point of the present 




Figure 1: Feynman diagrams calculated in this work. 

the approximation scheme: From eqn. (^) follows, that the identities ( p9|) should 
hold in every order of the resummed perturbation expansion in powers of S^. 

We now have to determine the three unknown parameters '^Fi'^t and 75;. For 
that purpose we calculate the self energy of the quarks and the polarization tensor 
for the gluons to one- loop order with the effective propagators. As is known from 
the exact integral equations for Green's functions, we therefore neglect only three- 
point and higher correlations. Up to this level however, our approximation scheme 
is completely consistent. 

As the parameters 7f,7t and 7^; are due to interactions they are necessarily 
proportional to some positive power of the coupling constant g times the temperature 
T. At temperatures high enough in order to ensure a small coupling constant the 
width parameters are necessarily small compared to the temperature. Therefore, 
we consider the limit of high temperature (for details see appendix B). As we want 
to provide a well-defined infrared limit with our parameterization of the spectral 
functions we determine 7f,7t and 7^; precisely in this limit. 



equation 



If = tTi [7oS^(po,0 



(31) 



The expression for the one-loop self energy S^(po,p) is, according to diagrammatic 
rules for the bottom diagram in figure |l|, 



E (p) = ig 



(32) 



For high temperatures and for positive we find, using the procedure outlined in 
appendix B: 

g^T^ Nl - 1 



1 



OiT) 



Nc (Po + 7t + If) 
For the electric and magnetic gluons the equations determining 'Je and jt are 



^71= noo(po,0)| 



Po^O 



and 



7t 



1 

2' 



2 - >J,n^,(po,o) 



(33) 



(34) 



(35) 



The gluon polarization tensor has a rather complicated structure. We therefore 
show only the electric and transverse part which we need for ( ^4]) and (pSf ); the 
following integrals correspond to the gluon tadpole, gluon/ghost loop and fermion 
loop diagrams depicted in figure |l|: 



IIe{po,0) = noo(po,0) 



(2^)4 
(2^)4 



^\2ko-fE + ll)SE{ko, k)SE{po + ko, k) 
+ 4(2A;o7t + 7t + k^)ST{ko, k)ST{po + ko, k) 
SFfi{ko, k)SF,o{po + ko, k) 
- SF,i{ko, k)SF,i{po + ko, k) 



(36) 



nT(po,o) 



ipj;n^,(po,o) 



2 



cos^(6') ^ ,, 3-cos^(6')^ , ■ 
^ ^SE{ko, k) + -^±ST{ko, k) 



v?. — Art„h Arr^^ 



1 - C0S2(^) 



■SE{ko, k)SE{po + h, k) 



2 o 2xl + COs2(0) 



+(27t(A;o + Po) + 7t - 2A;oPo - 3po 
■SE{ko,k)ST{po + ko, k) 



+ Ak' 



,1-COS2(0) 



S't(A;o, k)STipo + ko, k) 



Spfiiko, k)SFfl{po + ko, k) 

+ SF,i{ko, k) Pi Spjipo + ko,k) pj/p^ 



(37) 



Here we have used the simple relations between Se{p), Sl{p) and Sgh{p) due to 
The resulting expressions for high temperature are: 



nE(po,o) 



g^T^N, , g^T^Nf^F , g^T^N,{Q^T + Po) 



6 



+ 



nT(po,o) =po 



3(27,. + po) 
g^T^ ( iV, 



+ 



6(27t + po) 



+ 0(T) 



+ 



+ 0{T) . 



(38) 



(39) 



^27ir + po 27T+P0, 

Note, that quite similar results are obtained in a semiclassical approximation to 
the bremsstrahlung problem in hot and dense matter [^]. This indicates, that the 
high temperature limit we consider is closely related to a classical self-consistent 
calculation. 

Obviously, 11^(^0, 0) vanishes for po — suggesting 7^ = 0! On the other hand, 
it starts linear with po supporting a term 2po7T as in our ansatz for the transverse 
gluon propagator. We interpret this behavior as an inconsistency of our ansatz for 
the transverse propagator requiring a modification. 

The simplest modification of our ansatz such that St{po, 0)~^ vanishes for po ^ 



can be achieved by a spectral function like 
modified transverse propagator then reads: 

St{p) 



but now with Up 



p2 - ^2 _ rj.^^ 



(40) 



boP + 2|po|7t + P^ ' 

It seems that this propagator still has an infrared problem. However, recall that 
the proper quantity to consider are correlations of the field strength tensor rather 
than the vector potentials. In particular, the time averaged screening of an external 
magnetic field is given by the momentum space correlation function 



dt{Btit,x)B]{0,y)) ^ jd'pexp{ipx)Vf^p'STiO,p) cx Six) . 



(41) 



Here, we have not displayed any factors which are irrelevant for the spatial behavior. 
Obviously, static magnetic fields are screened and, in particular, the screening length 







iVc = 3 


Nf = l 


Nf = 2 


Nf = 3 


iV^ = l 


Nf = 2 


Nf = 3 


^^llE/{gT) 


0.913 


1 


1.080 


1.080 


1.155 


1.225 


It/ (gT) 


0.273 


0.323 


0.390 


0.315 


0.344 


0.380 


If/ (gT) 


0.199 


0.185 


0.168 


0.28 


0.271 


0.260 



Table 1: Self-consistent spectral width parameters for high temperature QCD 

Doing the previous calculations with the new spectral function for the transverse 
gluons yields the same results for the leading terms. A possible explanation for this 
feature is that at high temperature the kinetic terms are the dominant ones; these 
are unchanged by our modification. As the results are the same, our ansatz for the 
transverse gluons now is compatible with the infrared limit of the corresponding 
transverse polarization. The new equation determining 7^ reads: 

7T= f7^Pj.n^.(po,0)l . (42) 



4pi 







Po^O 



Now we can solve the set of equations (|3TD,(|3^) and (|42|) using the expressions 



(33), (38) and (B^). The solutions with positive values for the spectral width param- 



eters are 



7e = gTi 



lNf + 2Nc 



1 N^-l {gTf 
16 Nc 7f 



1 / 9 



IF = gT-^llN,- — -Nj-^/iNf-2N,)^ + 36iN^-l). (43) 



The solution for 'Je seems to indicate a gauge dependence. However, only for po = 
the electric propagator Se{p) is gauge invariant. There it reads 

1 

+ ^1% 

Obviously, ^7!; is independent of ^ and is just the famous "Debye" -screening mass 
for electric gluons 0, |22[] . 

Another remarkable feature of these results is the decrease of 7i7' with Nf. Indeed, 

for 

A^^ _ 1 







Nf = l 


Nf = 2 


Nf = 3 


ujT/{gT) 


0.457- 0.255 i 


0.496 - 0.263 i 


0.524 - 0.258 i 


ujtI (gT) 


-0.434 i 


-0.489 i 


-0.6 i 




iVc = 3 


Nf = l 


Nj = 2 


Nf = 3 


ujt/ {gT) 


0.541 - 0.309 i 


0.579 - 0.323 i 


0.614 - 0.331 i 


iVr/igT) 


-0.571 i 


-0.584 i 


-0.617 i 



Table 2: Transverse plasma frequencies of high temperature QCD 

we get a complex 'jp. This is not a very serious problem for the values of Nf and 
that are realized in nature; for A^^c = 3 the limit is Nf < 12. However, this behavior 
has to be understood from a purely theoretical point of view as it signals some 
strange and unexpected behavior. A possible explanation is that for large Nf the 
one-loop diagrams are not the appropriate ones to describe the leading behavior. 
In this case, one would have to consider a large-A^/ expansion instead of a loop 
expansion. 

In table we provide a listing of jEylT and from (^) in units of gT for 
several values of Nc and Nf. 

For completeness we also mention the results for pure Yang-Mills theory: 

It = gT—^ . (46) 
6 

One may also determine the plasma frequency of the system. We concentrate on the 
transverse oscillations, because only these are gauge invariant. They are solutions 
of the equation 

[iur)^ = UT{iujT,0) . (47) 

There is always one solution with up = because 11^(^0,0) is proportional to po- 
Numerical values for the other solutions are listed in table 0. A remarkable feature 
is the strong damping of these modes. Their imaginary part is about half of the 
real part. Moreover, there is also an overdamped mode. The ratio between the 



2 in contrast to the standard value -\/3 The explanation for this discrepancy 

may be the inclusion of an imaginary part of ut which is not negligible. 



4 The partition function of hot QED 



Essentially the same procedure can be applied to QED. Let us briefly recapitulate 
the most important steps. The partition function is now 



Zqed = j VA^VtpVipVcDc bxp^-Sqed) 



(4J 



with the action 

Sqed 



d X 



i'lixh^, {d^ - ieA^{x)) ipiix) + 



+-{K^A^{x)f + K^c{x)df,c{x) 



where / is again a "flavor" index. The electromagnetic field strength tensor is 



pa 



dfxA*^ — dyA"^ , 



(49) 



(50) 



and the gauge is the same as for QCD. 

We now define a modified unperturbed action Sq: 



^QED 



So 



So + Sint 

So + Sj + Sint ~ Sj 
So + Sint 
So + Sj . 



(51) 



Repeating the same procedure as for QCD means that we also introduce a spectral 
width for the ghosts. One might ask, whether this is not a contradiction to the 
fact that the ghosts are free fields that decouple from the physical fields. However, 
this question is irrelevant in our method, since we are dealing with generalized free 
fields having a continuous mass spectrum. Hence, we may choose any spectral width 
parameter for unobservable fields 

Our special choice for the ghost propagator is obtained by enforcing the BRST 
identities (^). This ensures the (gauge dependent) choice (p) = 0, and therefore 
is motivated by simplicity of computation. 



The equations for 'Jfj'Je and 7t are the same as in the Q CD-case, namely 



and P^). The fermion self energy now is 



E^(p) =ie^ 



Spip + k){2ST{k) + SL{k) - SE{k)) 



(52) 



For high temperature and for positive po "we find 

^Tr [roS^(po,0)] = ''f , ^ \ , + (^(r) • 

4 ' '8 (po + 7t + 7f) 



(53) 



The electric and transverse component of the photon polarization tensor read: 

nE(po,o) = noo(po,o) 

f d^k \ 

hi SFfi{kQ,k)SFfl{po + ko,k) 

- Sp^ko, k)SF,i{Po + h, k) 



-Ae^Nf 



(54) 



^t{Po,0) 



-pJ,n^,(po,o) 

o f d k 
ie^Nf 



(2vr) 



SF,o{ko, k)SFfi{po + ko, k) 



+ SF,i{ko, k)piSFj{po + ko, k)pj/p^ 
The resulting expressions for high temperature are: 

nMp„.o)^|;f^ + o(r). 

3(27f + po) 
UFipo,0)^,^^.OiT). 



(55) 



(56) 



(57) 



9(27f+Po) 

Solving the self-consistency equations for the spectral width parameters in QED 
then yields 



7b 



7t 



-'1 



{eTfNj 



367/ 



7f = ^\h;-^ 



eT 1 
YV2 



9 



(58) 



Again we find that for Nf > 4 the fermionic spectral width 7^7 becomes imaginary 



Finally, we also give the result for the transverse plasmon mode: 



1/5 1 

'^T = eT-^-Nf---i-fF (59) 
For standard electrodynamics with Nf = 1 this yields 

ut = eT^^(l - iv^) ^ eT (0.118 - 0.312i) (60) 
2 V 18 

This is a strongly damped oscillation. Here, the ratio between electric screening 
mass a/^7£ and the real part of ujt deviates even more strongly from than for 
QCD. However, since the plasma oscillation is also stronger dameped than in the 
QCD case, this only supports our previously stated conclusion that one may not 
negelect the imaginary part for the determination of the plasmon mode. 



5 Summary and conclusions 

In the present paper, we have investigated the self-consistent determination of finite 
temperature spectral broadening for QCD and QED in the infrared limit, i.e., for 
zero "particle" momenta. For the "electric" part of the gauge boson propagator, we 
reproduce the well-known Debye screening mass, whereas for the fermions as well as 
for the transverse gauge bosons we acquire new results: In contrast to the method 
of hard thermal loops we obtain, that all spectral width parameters •jp, 1e, It 
are of the order gT, respective eT, in the high temperature limit. 

Based on results for massive fermions we believe that in some regimes the 
constant spectral width obtained here might be a good approximation to the full 
solution of the problem. However, we have formulated a consistent recipe for the 
perturbative expansion of QCD and QED in terms of generalized free fields ||19|| , 
whose continuous mass spectrum is described by our effective propagators (eqn. 
(|30D, transverse boson propagator replaced by (^OD ). 

The skeleton expansion using these propagators may serve to calculate the full 
QCD and QED Green's functions [^, |]. Indeed, if the estimate (|^) holds for any 
infrared regulator, i.e. if M may be replaced by a typical 7, our results suggest that 
now a perturbative calculation of high temperature QCD or QED is possible; the 
expansion parameter in this cased would be g instead of g^. 

The virtues of a finite temperature perturbation theory using such propagators 
with a nontrivial spectral function lie at hand: 

1. It takes into account the Narnhofer-Thirring theorem described in the intro- 



2. It takes into account that the basis of a perturbative expansion has to be 
chosen according to the symmetry of the problem: The symmetry group of 
a temperature state is not the Poincare group with its stable particle rep- 
resentations, but rather possesses only representations with continuous mass 
spectrum |[. 

3. It is free of unphysical infrared divergences in every order, because products 
of mass-shell 5-functions cannot appear. 

Finally, our results for the transverse plasmon modes indicate a rather strong 
damping in contrast to the frequently used assumption that damping of these col- 
lective modes is suppressed by a factor of g or e. Due to the large imaginary part 
of these plasmon frequencies the relation rriei = v^Re(to'T) no longer holds. 
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A Notation 

We use a euclidean metric g^^, = 5^^. Space-time integrals are denoted by 



The euclidean 7-matrices are Hermitean and obey the anticommutation relations 



The gluon fields = A^T*^ are local elements of the Lie algebra SU(N). As such 
is traceless and Hermitean. The commutation relations of the generators T"" are 





(61) 



(62) 




(63) 



They are normalized such that 



Tr [T^T^] 




ab 



(64) 



and their completeness relation reads 




1 




(65) 



a/3 pa 



N 



This implies 

OT. = ^^^^- ■ (66) 
A useful relation for the structure constants reads 

Fourier transforms of a field (f){x) are defined as follows: 




Here, are the Matsubara frequencies for bosons (fermions) depending on the 
character of the field 0, and defined in (^. 



B High temperature limit 

Sums over Matsubara frequencies are transformed into contour integrals using the 
well-known identity: 

The upper sign is for bosonic z/^, the lower sign for fermionic z/„. 

The remaining one-loop integrals have the following typical form (see ref. ^ for 
details) : 

oo oo oo 

dE ldkeA,{E,k) JdE'^^^p^^ (70) 



exp ( — ) ± 1 







The E' integration as well as the /c-integration are carried out analytically, as contour 
integrals in the complex plane of the corresponding integration variable. 

For the remaining i?-integration, we use the fact that the integrand is an odd 
function of E depending only on the scale 7, 

f{E) = lf{-) , f{E) = -f{-E). (71) 

7 

multiplied by a temperature distribution function. For simplicity we discuss the 
case where f{E) has energy-dimension 1. 



The i?-integration may be split up into two parts, 

'',E {'f> = ,T L 4^ + 7^ k /';' . (72) 



exp(f)±l «^pW±1 / exp(4)±l 

In the first part, we may use an approximation for / at large arguments. 

The contribution of the leading term to the complete expression is of order T^. The 
remainder is at most of order 7T. 

In the second piece of the integrals, we expand the distribution function for small 
X. The resulting integrals are infrared finite because f{x) is finite at x ^ 0. As it 
is an odd function of its argument it vanishes at least proportional to x for small x. 
The leading term arises for bosons where the distribution function behaves as 

^ — + 0(1) . (74) 



exp (xj;^ - 1 

This yields a contribution of order 7T. 
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